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ABSTRACT 


Climate change has resulted in an increase in naval activity in the Arctic Ocean. While 
extensive research has been done on the interaction between different forms of ice and 
conventional steel hulls, little work has been done to assess the potential for composite¬ 
hulled vessels to operate safely in Arctic sea lanes. The goal of this experiment was to 
determine how well a composite hull can withstand ice loading and impacts with free- 
floating ice and whether a tumblehome, conventional, or a vertical hull performs best in 
the presence of ice. In order to do this, the tow tank at the Naval Postgraduate School was 
used. An ice equivalent was found which allowed for repeatable experimentation. Strain 
gauges were attached to critical locations of the composite plate towed through the tank. 
Both plate deformation and the resulting hydrodynamic force were measured for various 
ice conditions. A high-speed camera was used to qualitatively understand the flow 
characteristics of the ice around the composite plate. Finally, Ansys was used to 
determine if it was possible to replicate the experimental results. This work lays a basis 
for future tow ta nk experimentation involving ice interaction at the Naval Postgraduate 
School. 
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I. INTRODUCTION 


A, ARCTIC OPERATING ENVIRONMENT 

The emergence of sea lanes throughout the Arctic region is an on-going event that 
will have an effect on international politics for years to come. Not only is the region 
experiencing heavier maritime traffic due to more efficient shipping routes, but it also 
may become an area of interest for many nations in the surrounding area. Both the 
Beaufort and Chukchi seas are believed to hold the world’s largest unrecovered reserves 
of oil and natural gas [1]. In 2012, the U.S. Coast Guard began operation Arctic Shield in 
which it assessed its capabilities of operating in the aforementioned seas as Shell prepares 
to tap into their natural resources [1]. The U.S. Navy also conducts ICEX exercises, 
which familiarize submarines with operational procedures for transiting the Arctic. These 
are but a few of the actions the U.S. has begun to prepare for Arctic operations. While the 
initial probing of the natural resource rich seas in the Arctic has been peaceful, future 
global demands of such natural resources could spur conflict. 

Such a conflict may be similar to the current dispute of the Spratly Islands and the 
surrounding areas in the South China Sea. There have been numerous incidents involving 
vessels from China, Taiwan, the Philippines, and Vietnam as they each seek to lay claim 
over the natural resources in the region [2]. The problem lies in the United Nations 
Convention on the Law of the Sea, or UNCLOS. This signed document defines a 200 
mile exclusive economic zone (LLZ) that extends from a countries shoreline. The 
proximity of the countries in the South China Sea, shown in Ligure 1, skews each 
country’s claim to the islands. 


1 




CHINA 


TAIWAN 


China's claimed 
territorial waters 


UNCLOS 200 nautical mile 
Exclusive Economic Zone HAINAN South 


Disputed Islands 


China Sea 


Paracels 


VIETNAM 

f 

t 

I 

t 


Scarborouqh 
Shoal , 


Spratlys 


/ PHILIPPINES 


MALAYSIA 


,«RUNEI 

0 

MALAYSIA 


INDONESIA 


SOURCE: UNCLOS jod CIA 


Figure 1. South China Sea Geography Resulting in Territorial Disputes, 

from [2], 


In the 1970s and 1980s, incidents between Vietnam and China resulted in the 
death of over 130 Vietnamese sailors [2], Over 40 years later, the dispute over natural 
resources continues as each country attempts to validate its claim over the islands. 

While there have been no known clashes over territory in the Arctic, the same 
potential for disputes as found in the South China Sea is present. Canada, Denmark, and 
Russia all assert that the Lomonosov Ridge is an extension of their continental shelves 
[3]. The clustered Arctic region and various countries’ territorial claims are displayed in 
Figure 2. 
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Figure 2. Arctic Region Geography Showing Potential for Disputes over 

Natural Resources, from [4]. 


For this reason, it is in the United States’ best interests that it should be well 
versed in operating in the Arctic. The U.S. military has taken initial steps to prepare for 
Arctic operations. The DOD has published Arctic Strategy and the U.S. Navy Arctic 
Roadmap [3]. These documents promote the extension of infrastructure and additional 
research to support Arctic operations. ONR funds facilities such as NRL to conduct 
operations which better prepare the U.S. for Arctic operations. Countries such as Russia, 
Canada, and Norway have already begun honing their military prowess in the region [3]. 
This paper intends to add to the knowledge base by addressing the interaction between a 
fast-moving composite hulled ship and free-floating ice found in Arctic sea lanes. 
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B. 


COMPOSITE STRUCTURES 


Considering extensive research has been conducted on the interaction between 
traditional steel hulls and ice [5], this paper focuses on the interaction between E-Glass 
and ice. Composite material is becoming more prevalent in construction of cutting 
edge naval vessels. In 2006, the U.S. navy commissioned an experimental vessel named 
M80 Stiletto [6]. This was the navy’s first ship constructed entirely of composite 
materials. The 24-meter long ship is built from carbon fiber reinforced epoxy. The ship is 
displayed in Figure 3. 



Figure 3. Vessels like the Fully Composite Constructed M80 Stiletto have 
Potential to Operate in the Arctic Region, from [7]. 


The ship is intended for littoral operations. The carbon fiber hull results in a 
lightweight craft with a draft of one meter at its top speed. It is capable of speeds of up to 
60 knots [6]. As the craft travels faster, the hull is lifted farther out of the water, reducing 
the drag which inhibits forward motion. In conjunction with its shape, the composite hull 
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allows the M80 Stiletto to operate undetected to radar. These attributes make the M80 
Stiletto ideal for any special operation mission. 

Another topic of interest covered in this thesis is the performance of different hull 
types as they travel thorough the Arctic. More specifically, these hull types’ ability to 
travel through ice will be assessed. A hull is characterized by its profile at the waterline. 
The three hull types under investigation are the flared conventional hull type found on the 
Arleigh Burke Destroyer, the swept tumblehome hull found on the DDG-1000 Zumwalt 
Destroyer, and the wall-sided (termed “vertical” for the experiment) hull found in some 
tankers and transport vessels. These hull types are displayed in Figure 4, Figure 5, and 
Figure 6, respectively. 



Figure 4. Tumblehome Hull Exhibited by the DDG-1000 Zumwalt Destroyer, 

from [8]. 
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Figure 5. Conventional Hull Exhibited by the DDG-51 Arleigh Burke 

Destroyer, from [9]. 



Figure 6. Vertieal Hull Exhibited by a Medium Sized Norwegian Cargo 

Vessel, from [10]. 
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Most nautical vessels exhibit the conventional hull of the Arleigh Burke 
Destroyer-shaped hull. These hulls are most notable for their overall eonsistent 
performance and their stability in rough seas; henee every warship in the U.S. Navy other 
than the DDGIOOO has this hull type. The main purpose of the tumblehome design is to 
signifieantly reduce the radar cross section of this ship. Another benefit of a tumblehome 
hull is that the ship travels through the ealm seas more efficiently, which helps reduee the 
cost of fuel. These advantages eome at the detriment to the ship’s ability to operate in 
rough seas. Most crities believe the DDGIOOO cannot operate safely in hazardous 
eonditions. In 1997, NFS student Brian Ellis concluded that a tumblehome hull has the 
potential to reduee the maximum righting arm of a ship by 40% when eompared to 
conventional hulls [11]. The eurrent experiment will assess the ability of each hull type to 
displace and effectively travel through free-floating ice. 

C. PRIOR RESEARCH 

This study focuses on the impact between floating ice and a eomposite hull. 
However, as mentioned before, only one commissioned U.S. Navy vessel utilizes a 
composite hull. Composite materials are typically only found in the superstructure of 
large sea-going vessels. The bulk of researeh involving the impact resistance of 
composite materials stems from the aviation industry. Pinnell and Sjoblom [12] state that 
the aviation industry is eoneerned with two different types of impaets; high-velocity 
impacts which result in penetration of the composite and low-velocity impacts resulting 
in visual damage to the composite strueture. Their paper focuses on the results of low- 
veloeity impacts which are similar to those experienced by a ship’s hull as the hull strikes 
floating ice. Pinnell and Sjoblom found that low-velocity impacts can result in matrix 
cracking, delamination, and fiber fracture in two different loeations. The two loeations 
they described are the immediate layers at the location of the impact and in the plies just 
behind the baek surfaee of the composite plate [12]. A NASA teehnieal paper [13] 
supports these findings and notes that the most severe eracking and delamination oecurs 
on the side opposite of the impact point where the flexural stresses are greatest. Both 
visual inspection and X-Ray imaging were used to loeate surfaee and subsurface failures 
in the composite [13]. 
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A common method of testing the impact resistance of eomposite materials is by 
eompleting a drop tower test. This involves the release of a mass from a known height 
above a composite plate. This test produees repeatable test results. Russell’s [14] work 
shows the repeatability of subsequent trials using a drop tower. Russell’s work shows that 
the same impact parameters results in a different response for a composite plate in air and 
a composite plate in water. The maximum strain occurs in a different location of a plate 
and at a higher magnitude in water. This is due to the fluid structure interaetion (FSI). 
The density of water is nearly 784 times denser than air. This results in structures having 
different responses when perturbed in the different fluids [14]. Another work shows that 
the presence of water around the structure results in an inereased stress response from the 
impact of nearly 50% [15]. These results will be taken into consideration for this tow tank 
experiment. Ideally, a eomposite hull would be fabricated in whieh only one side would 
be exposed to water. However, a simple eomposite plate was used with water present on 
both sides of the plate as an initial study. 

D, OBJECTIVES 

The goal of this thesis was to eharacterize the response of various hull types as 
they impact free-floating iee which is commonly found in Arctic sea lanes. In order to 
perform repeatable tow tank experiments, an ice equivalent was found. Upon proof of 
coneept, this ice equivalent may be used in future tow tank experiments which hope to 
expand the knowledge base for Arctie operations. Video of the front edge of the 
eomposite plate was recorded to determine whether the ice equivalent strikes the plate or 
is displaeed by the bow wave generated by the plate. It was possible that there was a 
eritieal veloeity at which the ice is no longer displaeed by the bow wave. Force and strain 
measurements were reeorded to quantify how well each hull orientation progressed 
though iee infested waters. Two different plates were used in testing. The two plates were 
same size; however, the boundary conditions of the two were different. The first plate 
was ereated for researeh in a prior thesis. This plate has a frame surrounding the plate. 
This simulates Diriehlet boundary eonditions around the plate. Another eomposite plate 
was fabrieated and used to simulate Neumann boundary conditions. 
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A computer model of the tow tank experiment was also generated in hopes of 
replieating the experimental data. A two fluid domain consisting of air and water was 
ereated in CFX to assess the free surfaee effeet at the water’s surfaee and its interaetion 
with the three different hull orientations. A base model was generated which could be 
improved to better simulate the interaetion between free-floating objects and the 
composite plate. 


9 



THIS PAGE INTENTIONALLY LEET BLANK 


10 



II. EXPERIMENT 


A. TOW TANK 

This experiment utilized the tow ta nk located in Halligan Hall on the Naval 
Postgraduate School campus. The tow ta nk and its dimensions are shown in Figure 7. 



Figure 7. Halligan Hall Tow Tank with Dimensions. 


In order to tow objects, a carriage is attached to the top rails with ball bearings. 
This allows the carriage to translate along the length of the tow tank. The carriage is 
pulled using a five horse-power motor connected to a pulley system which runs along the 
length of the tow tank. This motor is controlled with a Baldor controller connected to 
30 meters of cable. This allows for the operator to control the carriage at any location 
around the tank. The operator inputs speeds settings in Hz. A minimum of 2.3 Hz is 
required to move the carriage. The full set of operating guidelines can be found in the tow 
tank SOP located in Appendix A. The relationship between the input frequency and the 
carriage velocity will be shown later. The controller is shown in Figure 8. 
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Figure 8. Tow Carriage Controller. 


In order to tow the composite plate at a specified angle through the water, an 
angle selector is located on the underside of the carriage. This part was fabricated by the 
machine shop for previous theses. The angle selector is shown in Figure 9. 



Figure 9. Angle Selector for Precise Plate Orientation. 
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There are two bolt holes in the attachment bar, which align with the top hole and 
the bottom arc of the angle selector. The attachment bar is a rectangular aluminum tube 
that connects the towed composite plate to the tow carriage. After tightening the top bolt, 
the bar can be moved along the arc until the intended plate angle is obtained. The bolt 
along the arc is then tightened to lock the plate orientation in place. This experiment was 
done for a vertical plate and with the plate at ±30°. 

The pulley system is connected to the front and rear of the carriage. The pulley 
wire is bolted to the rear of the carriage in a solid connection. Located on the front of the 
carriage is a Honeywell Model 41 “Pancake” type load cell. This load cell measures the 
force at which the pulley system tows the carriage. The other data obtained from this 
experiment was the strain experienced by the composite plate. The strain gauges were 
connected to quarter-inch bridge adapters. These adapters and the front load cell were 
connected to a National Instruments Wireless Sensor Network (WSN). The WSN only 
allowed for four channels of experimental data to be streamed for data storage to a laptop 
using a Wi-Fi connection. The four channels were occupied by three strain gauges and 
the front load sensor. The procedure for connecting the strain gauges to the quarter-inch 
bridge adapter can be found in Appendix B. The necessary equipment for data acquisition 
shown in Figure 10 can be found fixed to the top of the tow carriage. 
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Figure 10. WSN (A), Quarter-Inch Bridge Adapters (B), and External Battery 
Pack (C) Used to Transmit Strain and Force Values for Data Storage. 


The experimental data was recorded to the laptop using a program called 
Labview. The strain gauge factor was entered into the program in order to obtain valid 
strain measurements. The computer program display is shown in Figure 11. Example 
force data from an experimental trial is still displayed on the strip chart. 



Figure 11. Fabview Computer Program Used for Data Acquisition. 
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The gauge faetor for each strain gauge is entered into the appropriate “GF Ch #” 
for each strain gauge. The default sample rates were used at 100 samples per second. 
Real time data can be viewed in the top strip chart. After each trial this chart was cleared 
of data. Selecting the “Save to File” button allowed for the experimental data to be saved 
as a tab delineated fde which could be imported to excel for data reduction. 

The data for each trial can be divided into three stages. The first is the 
acceleration stage. In this stage, the force peaks rapidly as the pulley system applied the 
initial force required to begin moving the carriage. The next stage is the steady state 
condition. At this point, the data is consistent and oscillates about an average value. This 
is the stage of interest. All data reduction will be conducted in this range of data. Last is 
the deceleration stage. This stage has the similar erratic behavior found in the 
acceleration of the carriage and will not be used for data reduction. The typical data for 
one trial is shown in Figure 12. 



Figure 12. Data Obtained from Experimental Trials Divided into Acceleration 
Stage (A), Steady State Stage (B), and Deceleration Stage (C). 


The data has been divided into the three stages. For each run, the steady state 
range will be visually selected to match the steady state behavior displayed in Figure 12. 
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The final data required for the experiment was the veloeity of the earriage as it 
towed the eomposite plate through the water. This was obtained using an Olympus 
i-SPEED high-speed eamera. In addition to the veloeity measurement, the eamera 
allowed for the flow characteristics of the ice equivalent around the plate to be assessed. 
The camera was placed in front of the glass viewing pane of interest. The average 
velocity was determined by dividing the distance the plate traveled in the camera’s field 
of view by the total elapsed time required to cover the distance. 

B, ICE EQUIVALENT 

Ideally, this experiment would be done in an ice tank. An ice tank is a 
temperature-controlled facility that houses a towing tank. The temperature control allows 
the working area to remain below the freezing temperature of water, allowing the use of 
real ice for experiments. There are currently 11 ice tanks around the world. The largest is 
102 meters long, 10 meters wide, and 2 meters deep which can be found in St. 
Petersburg, Russia. Since NPS does not have a climate controlled area, using ice was not 
possible. Eor this reason, a viable substitute had to be found which allowed for 
repeatability in all of the experimental trials. The main goal of this experiment was to 
match the density ratio between the fresh water tow tank water and the simulated ice to 
Arctic brash ice and Arctic Ocean seawater. This would give the most realistic floating 
characteristics for the ice equivalent for the experiment. Brash ice is the accumulations of 
floating ice made up of fragments not more than 2 meters across [16]. An example of a 
ship navigating through brash ice is shown in Eigure 13. 
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Figure 13. Vessel Transiting through Brash lee, from [17]. 


1, Research 

To properly simulate Aretie iee, it is first neeessary to have an understanding of 
the Aretie environment and the interaetion between the iee and sea water. When sea 
water freezes, the salt and the water separate into two different phases: brine and iee. If 
the salt is eaptured inside the body of iee, it remains as a poeket of brine. Brine is a salt 
water mixture with a very high eoneentration of salt. Any salt that is not physieally 
eontained inside the iee dissolves back into the Arctic Ocean [18]. This characteristic 
results in sea ice having a large range of density values. Numerous density measurements 
of sea ice have been compiled by Frederking and Timco [19] and are shown in Table 1. 
The density of sea ice ranges from 720 kg/m to 960 kg/m . Most recorded values are 
nearest to the upper limit of that range. 
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Table 1. Density Measurements of Sea lee, from [19]. 


Invesligaior 

Measurcmcni 

technique 

Density range 
Mgm ’ 

Density 

(above waterline) 
Mg m ' 

Deiksit) 

< helow waterline) 
Mg m ’ 

Brine 

drainage 

Fmt\far sea ict 

Maliiigrcn (1927) 

- 

0.91410 0.924 

- 

091410 0.924 

NO 

Butkiivich (1956) 

mass/volume 

0.78 to 0.92 

0.85 to 0.92 

0.88 to 0.92 

YES 

Weeks and Lee (I9.')8) 

mass/volume 

0 77 toO.% 

- 

0.86 to 0.91 

YES 

■.angkben (19,59) 

mass/volume 

0.84 to 0.90 

0.84 to 0.90 

090 to 093 

NO 

Pixmder and l.iille (I9S9) 

specihe gravity 

0.81 to 091 

- 

- 

YFS 

Pounder and Sulinski (1960) 

spcciHc gravity 

0.943*0.02 

-- 

- 

YES 

Brovin (1962) 

specific gravky 

092 

- 

0.92 

NO 

Ijuiglebcn and Pounder (1962) 

displacemeni 

0.884 to 0.91 

0.884 

0.91 

YF.S 

Kuhnen (1972) 

mass/volume 

0.90 to 0.92 


0.90 to 0,92 

NO 

Nakawol 198.7) 

displacemeni 

0.88 to 0.92 

0.88 to 0.91 

0.90 to 0.92 

NO 

Timeo and Frederking (198,1) 

mass/volume 

0.9(M 

- 

0904 

NO 

Sinha (1984) 

mass/volume 

090 to 0.91 

0.90 to 0.91 

0 90 to 0.91 

NO 

Sinha (1986) 

mass/volume 

0.90 to 0.91 

0.90 U> 091 

0 90 to 0 91 

NO 

Timeo and Frederking (1986) 

mass/volume 

0.76 to 0.89 

- 

- 

YES 

Urube and Inoue (1986) 

mass/volume 

0.75 to 0.88 

0.75 

088 

YES 

MuUt sear rrr 

Malmgren (1927) 

- 

0.857 to 0.893 

0.857 

0,885 K> 0.893 

YES 

Langteben and Pounder (1962 ) 

displacement 

0.91 

- 

- 

- 

Cherepanov { 1966) 

- 

0 912 to 0.917 

0912to0915 

0.915 to 0.917 

NO 

Kovacs and .Mellor ( 1971 ) 

- 

0 84 to 0.94 

0.84 Ui 0.88 

0.90 to 0.94 

NO 

Hibler «i al. < 1972) 

- 

0 910 to 0.915 

- 

- 

- 

Sinha (19X4) 

mass/volume 

0.72 to 0.90 

0.72 to 0.85 

0,K9io0.90 

NO 

Jeffriesct al (1988) 

mass/volume 

0.87 to 0.912 

- 

- 

- 


With the density of sea lee known, it was then neeessary to find the density of 
Aretie Oeean seawater. The density of salt water is a funetion of its temperature and 
salinity aeeording to the Equation of State of Seawater [20]. In order to obtain values for 
temperature and salinity, the NOAA database was aeeessed. The database requires a 
season and a depth in order to display the expeeted data for salinity and water 
temperature. Sinee Aretie sea lanes would be most aeeessible in Spring and Summer, 
those times were seleeted. The surfaee layer was seleeted for depth. The results for the 
salinity and temperature are found in Figure 14 and Figure 15, respeetively. 
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Figure 14. Expected Surface Layer Practical Salinity Unit (PSU) in the Arctic 

Region Published by NOAA, after [21], 



Wi>rt<Ociaaw *aM n »fcw<l C* p>iloiDg » Aictc 
C u ^ lagi*.! 



Figure 15. Expected Surface Layer Temperature (°C) in the Arctic Region 

Published by NOAA, after [21]. 


This data shows that the salinity varies from 29-35 PSU and the temperature 
varies from -2°C to 6°C. The maximum extreme of the density value results from the low 
temperature, high salinity water. The opposite results in the minimum extreme density 
value. Using the Equation of State of Seawater, the density of the Arctic seawater values 
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ranged from 1022.8 kg/m to 1028.2 kg/m . This results in a final density ratio of sea iee 
to Aretie Oeean seawater ranging between 0.70 and 0.94. 

In order to get the eorreet floating behavior of the simulated iee in the fresh water 
tow tank, the density ratio between the iee and oeean water must mateh the simulated iee 
and the tow tank water. It was first assumed that the water used in the tow tank had the 
typieal density of 1000 kg/m eornmon to freshwater. This meant the simulated iee must 
have a density ranging from 700 kg/m to 940 kg/m . The material must also not absorb 
the water, as that would ehange its density during testing. After searehing various 
materials, two seemed best for the experiment. The first was High Density Foam (HDF) 
maehining board, whieh is eornmon in model making. The next was Low Density 
Polyethylene (LDPE) plastie. After researehing eost, resistanee to water absorption, and 
handiness, LDPE was ehosen as the material to simulate iee. LDPE eould be ordered in 
the shape of a rod whereas the HDE would have to be maehined into an iee shape from a 
large reetangular bloek. The rods simply had to be eut every few inehes to ereate the 
simulated iee. MeMaster Carr advertised a density of 913 kg/m for the EDPE [22]. This 
fit well in the desired range of ideal densities. 

2. Fabrication 

In order to fulfill the iee requirements for the desired experiments, two different- 
sized rods were ordered. Two-ineh diameter rods were ordered whieh would simulate a 
ship moving through a field of small randomly dispersed sea iee masses. Earger four-ineh 
diameter rods were purehased to simulate impaets with large bodies of iee. The EDPE 
rods are shown in Eigure 16. 
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Figure 16. LDPE Rods Used to Simulate Aretie lee. 


The density of the LDPE and the tow tank water was then verified to ensure the 
experiment properly represented the Aretie environment. Eirst, each rod was weighed 
using a large electric scale, as shown in Eigure 17. 



Eigure 17. Electric Balance Used in Density Measurements. 
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Then, the diameter of each rod was measured at the top, middle and bottom using 
calipers. The rod diameter is shown in Table 2. 


Table 2. Diameter Variations along LDPE Rods. 


Rod 

Number 

1 

2 

3 

4 

5 

6 

Top 

Diameter 

(mm) 

51.47 

51.49 

51.43 

51.53 

51.48 

51.44 

Mid 

Diameter 

(mm) 

51.46 

51.46 

51.43 

51.46 

51.43 

51.45 

Bottom 

Diameter 

(mm) 

51.43 

51.45 

51.44 

51.56 

51.39 

51.40 

Average 

(mm) 

51.453 

51.466 

51.433 

51.517 

51.433 

51.43 


Finally, the length of the rods was measured and their densities calculated. The 
final results are shown in Table 3. 


Table 3. Density Calculation Summary. 


Rod 

Number 

1 

2 

3 

4 

5 

6 

Mass (kg) 

2.3423 

2.3445 

2.3471 

2.3472 

2.35333 

2.35587 

Length 

(m) 

1.2192 

1.2200 

1.2208 

1.2224 

1.2256 

1.2287 

Average 

Diameter 

(m) 

0.051453 

0.051466 

0.051433 

0.051517 

0.051433 

0.05143 

Volume 

(m') 

0.002533 

0.002538 

0.002536 

0.002548 

0.002546 

0.002553 

Density 

(kg/m^) 

924.8 

923.8 

925.4 

921.2 

924.2 

922.9 


The final average density of the LDPE was 923.7 kg/m . This is slightly above the 
published density is 913 kg/m . A graduated cylinder and an electric balance were then 
used to determine the density of the water in the tow tank. This is shown in Figure 18. 
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Figure 18. Graduated Cylinder and Scale Used for Water Density 

Measurements. 


The tank was filled and allowed ample time to reach temperature equilibrium with 
the ambient air. Three different trials were done with water taken from the middle and 
outer edges of the tow tank. The density results are shown in Table 4. 


Table 4. Density Results for Tow Tank Water. 



Volume (ml) 

Mass (g) 

Density (kg/m ) 

Fill End 

93.0 

92.3 

992.4 

Middle 

98.8 

97.5 

986.8 

Drain End 

97.3 

96.3 

989.7 


Average = 989.7 


The final density ratio between the LDPE and the tow ta nk water was 0.933. The 
expected value of 0.913 for the density ratio between the LDPE and water would have 
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been more ideal for the experiment sinee it was eloser to the majority of reeorded sea ice 
density values; however the actual density ratio used in the experiment still fell into the 
range of density ratios actually recorded for sea ice. 

In order to have consistent floating characteristics, the rods were cut into various 
lengths and placed in the tow tank. The two-inch diameter rod was cut into segments 
ranging from one inch to two-and-a-half inches, at half-inch intervals. The resulting 
floating characteristics are shown in Figure 19. 


4 ' 



Figure 19. Floating Characteristics of Rods with Lengths 25.4 cm, 38.1 mm, 
50.8 mm, and 63.5 mm (1.0 in, 1.5 in, 2.0 in, 2.5 in). 


The 38.1 mm (1.5 in) long segment was the longest section that had self-righting 
properties. The rods were cut using a table mounted buzz saw. Since many cuts were 
required, a vertical piece of metal was C-clamped to the work bench at a distance from 
the buzz saw that ensured the rods were cut into 38.1 mm (1.5 in) segments. The setup is 
shown in Figure 20. This setup allowed the rod to be continuously fed into the buzz saw 
while also being cut to a near-perfect length. 
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Figure 20. Small Diameter LDPE Rod Cutting Setup to Maximize Efficiency. 

The larger four-inch diameter rods were cut using a vertical band saw. After 
completing a similar float test to the one mentioned above, the large diameter rod was cut 
to a length of 76.2 mm (3 in). 

C. SIMULATION OF OPERATIONAL ENVIRONMENT 

This thesis simulated both small-scale and large-scale collisions between a 
composite hull and free-floating ice. The small-scale ice was represented by 38.1 mm 
(1.5 in) long LDPE rod sections with a diameter of 50.8 mm (2 in). The large-scale ice 
was represented by 76.2 mm (3 in) long LDPE rod sections with a diameter of 101.6 mm 
(4 in). Multiple scenarios for each ice equivalent size were completed. 

1, Small-Scale Ice Interaction 

This experiment utilizes two different sizes of simulated ice in order to account 
for different ice conditions that a vessel may encounter in Arctic sea lanes. The first sea 
state involves a large quantity of small pieces of ice, as shown in Figure 21. 
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Figure 21. Small-Scale lee Present in Arctic Waters, from [23]. 

The two-inch diameter rods were used to simulate this environment. To ensure 
repeatability between experiments, a specific surface capacity of LDPE rods was used to 
determine how many rod sections should be placed in a section of the tow tank. The 
surface capacity was calculated by dividing the product of the number of rods and the 
area of each rod by the maximum area of the rectangular section in the tow tank that 
could be occupied by the circular ice equivalent. The length of the glass panels of the 
tank determined the length of the area occupied by the ice equivalent. Each glass panel is 
37 inches long. The glass panel allows for the use of the high speed camera to observe 
flow characteristics and the interaction between the towed plate and the ice equivalent. 
This experiment included low, medium and high ice capacity trials. The experimental 
setup is for the different surface densities are shown in Eigure 22. The green line 
represents the towed composite plate. 
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Figure 22. Experimental Schematic to Simulate Conditions in Figure 21. 


In order to accomplish this experimentally, the number of rods required to fulfill 
the surface density requirement were dumped into the tow tank. They were allowed to 
float randomly within the 37-inch long section of the tow tank. Enough time was given 
for the rods and tow ta nk water to reach an equilibrium state. 

2, Large-Scale Ice Interaction 

The second simulated Arctic operating environment accounted for large free- 
floating ice that may be present in sea lanes. The goal was to model the impact between 
one or two large pieces of ice with the composite hull. This operating area would be 
similar to that shown in Eigure 23. 
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Figure 23. Scattered Large-Scale Ice Present in Arctic Waters, from [24], 


Multiple possible impact situations were assessed using the large-diameter LDPE 
rods. These include a single centered strike, a single off-centered strike, and a double 
impact. The tow tank schematics in Figure 24 display the different impact situations. 
These scenarios were completed for the three different hull angles. Videos were taken to 
qualitatively assess the interaction with each hull shape and the LDPE rods. 
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Figure 24. Experimental Schematic to Simulate Conditions in Figure 23. 

D, COMPOSITE PLATE FABRICATION 

The composite plate used in the tow tank experiment was made using woven glass 
fabric, resin, and hardener. First, the woven glass fabric was cut to size. Since the plate 
would be subjected to notable hydrodynamic forces it was decided that the plate would be 
10 layers thick. The desired plate dimensions were 30.5 cm by 30.5 cm (12 in by 12 in). 
Since the woven glass has the potential to move slightly during the layup process, the 
woven glass fabric was cut into 10 squares with a side length of 13 inches using a sharp 
cutting wheel. The setup is shown in Figure 25. 
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Figure 25 . Cutting Wheel Used to Cut Woven E-Glass Sheet. 


Next, other necessary materials that were required for the composite plate layup 
were cut to size. These include release ply paper, nylon peel ply, perforated release ply, 
bleeder cloth, and a vacuum bag. With the exception of the vacuum bag, these were all 
cut into squares with a side length two inches longer than the woven glass fabric. This 
ensured that any excess epoxy would not adhere to the table on the bottom side and 
would be absorbed by the bleeder cloth on the top. With the necessary materials ready for 
layup, the resin and hardener were mixed according to the manufacturer’s instructions. 
Pro-Set M1002 Resin and Pro-Set 237 Hardener were used. The two were mixed for five 
minutes using a weight ratio of 100:24. These were mixed in a mixing container inside of 
a fume hood. The setup is shown in Figure 26. The table used for layup consists of a large 
glass pane resting on top of a work bench. The workbench used is shown in Figure 27. 
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Figure 26. Resin and Flardener Mixed by Weight Ratio in Mixing Container. 



Figure 27. Workbench Used for Composite Plate Layup. 

To begin, the glass work area was cleared using a metal scraper. The bottom layer 
consisted of one layer of release ply paper which prevented the plate from curing to the 
glass work table. Next, a base layer of epoxy was rolled evenly onto the ply paper using 
the metal roller shown in Figure 27. Alternating layers of glass fiber and epoxy were 
rolled until the ten layers were completed. Next, the nylon peel ply, perforated release 
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ply, and bleeder eloth were plaeed on top. The bleeder eloth absorbed any exeess epoxy 
applied during the rolling proeess. Finally, a large vacuum bag is sealed around the entire 
set up using double sided tape. The electrical pump was used to create a vacuum of 
15 mmHg. The plate was given 24 hours to cure. The plate was then cut to size in the 
machine shop. Holes were drilled where the plate was to be bolted to the tow carriage 
bar. 

E, APPLICATION OF STRAIN GAUGES 

Strain gauges were used to measure the strain experienced during each of the 
scenarios depicted in Section C. As mentioned before, the data acquisition system can 
only transmit data from three strain gauges. For this reason, an Ansys model was created 
to visualize the strain incurred by the pressure of the water against the front face of the 
composite plate. The strain gauges were placed in accordance with the results of the 
Ansys model. The strain gauges were then applied using a series of epoxy and protective 
coatings which allowed the strain gauges to operate in the presence of water. 

1. Determining Strain Gauge Location 

First, it was necessary to determine the critical location on the plate where the 
hydrodynamic forces cause the most deformation in the plate. The data acquisition can 
only record data from three strain gauges, so the location cannot be chosen haphazardly. 
In order to do this, a simplified static structural simulation was done in Ansys to obtain 
qualitative information regarding the strain profile resulting from moving the plate 
through water. The plate was created in Solidworks and imported into Ansys. The 
mechanical properties of the E-glass were input into the engineering data. The data was 
obtained from a prior NPS thesis that used the same E-glass [25]. A fixed support was 
placed on the back of the plate where the attachment bar is bolted in place. Eigure 28 
shows the simulation and the experimental setup. 
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Next, a pressure was added to the front of the plate at the intended water level for 
experiment. Since this simulation was done before testing, the added height of the bow 
wave preceding the plate was not accounted for. Figure 29 shows the relationship 
between the hydrodynamic pressure in the experiment and the ANSYS model. 



Figure 29. Hydrodynamic Pressure Resulting from Towing Plate. 


ANSYS was used to solve for the equivalent Von-Mises stress in the plate. The 
results are shown in Figure 30. In order to have a better understanding of the location of 
the fixed support, the setup image is shown next to the final results. 
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Figure 30. Location of Strain Gauges Noted by Red X’s on Stress Profile. 


This image shows the strain is greatest just below the attachment bar. Since wires 
had to be soldered to the strain gauges they could not be placed directly underneath the 
attachment bar. The top line indicates how far the attachment bar extended onto the plate, 
which was 11.18 cm (4.4 in) from the top of the plate. The bar was 2.54 cm (1 in) wide. 
The second line shows the plane on which the first strain gauge was placed. A strain 
gauge capable of measuring strain in the vertical direction was placed on the centerline of 
the plate at this location. One 90° Tee Rosette capable of measuring strain in both the 
vertical and horizontal directions was placed midway between the center and the side and 
midway from the first strain gauge to the bottom of the plate. The locations of the strain 
gauges on the framed composite test plate had an impact on the locations of the strain 
gauges as well. This allowed for the strain resulting from the two different boundary 
conditions to be assessed. 

2, Strain Gauge Application 

With the location of the strain gauges established, the plate was prepared for the 
bonding agent. The location of each strain gauge was just above the five-inch grid line. 
The procedure is listed below. 

1. Cleaning Strain Gauge Location; The two strain gauge locations were 
sanded with fine sand paper to ensure the bonding surface was clean of 
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any imperfections and any contaminates that might result in delamination 
of the strain gauge. The area was then cleaned with ethanol and a gauze 
pad. This cycle was repeated twice. 

2. Placement of Strain Gauges: Using tweezers, the strain gauges were 
moved into place. Scotch Tape was then placed over the strain gauges 
and then peeled back into a loop. This allows the bonding agent to be 
applied to the underside of the gauge while ensuring the strain gauge will 
remain in its intended location. The end result of this process is shown 
in Figure 31. 



3. Application of Bonding Agent: An M-Bond AE-10 Adhesive Kit 
manufactured by Micro-Measurements was used to adhere the strain 
gauge to the composite plate. The kit consists of a two part epoxy, a 
dropper, and a plastic stirring rod. The dropper was fdled up to the preset 
mark with hardener. The hardener was added to the small jar of resin. 
After finishing the required five minute mix time, a very small amount of 
epoxy was applied to the strain gauges. The tape was then unlooped and 
pressed against the plate firmly. 

4. Curing Time: Weights were placed on top of the strain gauges to ensure 
they remain flush against the surface. This is shown in Figure 32. While 
the epoxy does partially cure in six hours, it was given 48 hours to achieve 
its maximum elongation capabilities. This will ensure that when the plate 
is flexed from the hydrodynamic forces, the strain gauge will not 
experience delamination from the plate’s surface. 
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Figure 32. Application of Weights to Ensure Proper Bonding between Strain 

Gauges and Composite Plate. 


5. Soldering Leads; After the 48 hour cure time, the tape was removed. This 
exposed the terminals to which the wire that connects the gauge to the 
sensor is applied. After ensuring the terminals were clean with a dental 
pick, the leads were then soldered to the terminals. 

6. Waterproofing; There were two coats required to properly insulate the 
strain gauge and wire connection from the tow tank water. The first was a 
general purpose laboratory coating called M-Coat A. A liberal coating was 
applied over the entire strain gauge and exposed wire. This coat was 
allowed to cure for 6 days as recommended by the manufacturer. Next, 
Dow Corning MIL-A-461 RTV Coating, a room temperature curing 
silicone rubber, was generously applied on the strain gauge and wire 
connection. The wires were then routed through the attachment bar. This 
eliminates the need to tape the wires to the attachment bar. The 
waterproofed strain gauges are shown in Figure 33. 
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Figure 33. Waterproofed Strain Gauges. 


F. TEST PLAN 

The goal of this experiment was to characterize how well each hull performed in 
the presence of the different-sized ice equivalent. This test plan was completed with both 
the framed composite plate and the frameless plate. The framed composite plate is shown 
in Figure 34. The strain experienced by the framed plate is expected to decrease for all 
trials when compared to the strain experienced by the frameless plate. 



Figure 34. Framed Composite Plate Used in Testing. 

37 
















In order to compare performances of each hull type, the strain experienced by the 
plate and the force required to travel through the media were recorded for the base water 
trials, 50% ice capacity, and for an impact with two large piece of ice. Additional trials 
for smaller ice densities were completed to determine the relationship between increasing 
ice density and the resultant forces and strain values. The other situations described in the 
Ice Equivalent section were also tested. The full test plan followed to obtain the 
necessary data is shown in Table 5. Figure 35 shows the difference between the 
tumblehome orientation and the conventional orientation. As mentioned before, the two 
orientations will be angled at 30°. 


Table 5. Test Plan. 


Run 

Trial Type 

Desired Data 

Speed/Orientation 

1 



3 Hz / Vertical 

2 

Base Air Trials 

Force 

5 Hz / Vertical 

3 

7 Hz / Vertical 

4 



9 Hz / Vertical 

5 



3 Hz / Vertical 

6 



5 Hz / Vertical 

7 



7 Hz / Vertical 

8 



9 Hz / Vertical 

9 



3 Hz / Tumblehome 

10 

Base Water Trials 

Force, Strain, Velocity 

5 Hz / Tumblehome 

11 

7 Hz / Tumblehome 

12 



9 Hz / Tumblehome 

13 



3 Hz / Conventional 

14 



5 Hz / Conventional 

15 



7 Hz / Conventional 

16 



9 Hz / Conventional 

17 

16.67% Capacity Small Ice 
Trials 


9 Hz / Vertical 

18 

Force, Strain 

9 Hz / Tumblehome 

19 


9 Hz / Conventional 

20 

33.3% Capacity Small Ice 
Trials 


9 Hz / Vertical 

21 

Force, Strain 

9 Hz / Tumblehome 

22 


9 Hz / Conventional 
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Table 5 ( 

cont’d) Test Plan. 

Run 

Trial Type 

Desired Data 

Speed/Orientation 

23 

50% Capacity Small Ice Trials 

Force, Strain, Video 

3 Hz / Vertical 

24 

5 Hz / Vertical 

25 

7 Hz / Vertical 

26 

9 Hz / Vertical 

27 

3 Hz / Tumblehome 

28 

5 Hz / Tumblehome 

29 

7 Hz / Tumblehome 

30 

9 Hz / Tumblehome 

31 

3 Hz /Conventional 

32 

5 Hz / Conventional 

33 

7 Hz / Conventional 

34 

9 Hz / Conventional 

35 

Single Centered Large Ice 
Trial 

Force, Strain, Video 

9 Hz / Vertical 

36 

9 Hz / Tumblehome 

37 

9 Hz / Conventional 

38 

Single Off-Centered Large Ice 
Trial 

Force, Strain, Video 

9 Hz / Vertical 

39 

9 Hz / Tumblehome 

40 

9 Hz / Conventional 

41 

Double Centered Large Ice 
Trial 

Force, Strain, Video 

3 Hz / Vertical 

42 

5 Hz / Vertical 

43 

7 Hz / Vertical 

44 

9 Hz / Vertical 

45 

3 Hz / Tumblehome 

46 

5 Hz / Tumblehome 

47 

7 Hz / Tumblehome 

48 

9 Hz / Tumblehome 

49 

3 Hz /Conventional 

50 

5 Hz / Conventional 

51 

7 Hz / Conventional 

52 

9 Hz / Conventional 
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Figure 35. Conventional Hull (Left) and Tumblehome Hull (Right) with 0=30°. 

G. PROCEDURE 

Before eaeh trial, the earriage was positioned near the left end of the tow tank. 
This ensured that the carriage reached steady state conditions when it arrived at the center 
of the tank. This was important since the ice equivalent would be placed near the center 
of the tank. Using previous NPS tow tank experimental data [25], test speeds of 3 Hz, 5 
Hz, 7 Hz, and 9 Hz were chosen. According to previous data, these inputs should result in 
carriage velocities ranging from 0.40 m/s to 1.11 m/s. Approximately one m/s was 
desired as a maximum carriage speed for safety consideration due to the short length of 
the tow tank. The first step was to verify the speed results from the previous tow ta nk 
experiments. 

The speed of the carriage was determined using the high speed camera. Since the 
plate glass and the water refract the visible light in accordance to Snell’s Law, the 
distance the plate traveled in the video could not be perfectly calculated. For this reason, 
the distance traveled by the plate was determined experimentally. The camera was set up 
at the distance which allowed the entire glass pane to be in view. The plate was moved 
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just into view of the camera and stopped. The nearest vertical edge of the plate was the 
point of reference. The starting location was marked on the tow tank. Next, the plate was 
moved across the field of view of the camera and stopped at the far right of the camera’s 
field of view. At this point, only the rear bolts of the attachment bar were visible. The 
final location was marked on the tow tank and the distance between the two points were 
measured. Figure 36 shows the imaging process used to determine the distance traveled 
by the plate. 



Figure 36. Beginning and End Locations Used for Measuring Average Speed. 
Only the Attachment Bar Bolts are Visible as Indicated by the Arrow. 


It was determined that the carriage travelled 1.07 m (42.2 in) to pass through the 
camera’s field of view. The time stamp of the frame at the beginning and end of each trial 
was recorded to determine the average velocity of the carriage. The velocity of the 
carriage for each speed setting is shown in Figure 37. 
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•Speed (m/s) 
•Motor RPM 


Input Speed (Hz) 


Figure 37. Correlation between Input Speed Setting (FIz) and both Velocity 

(m/s) and Motor RPM. 


With the plate velocity known for the four desired speed inputs, the experiments 
described previously in the test plan were conducted. The base trials without the ice 
equivalent were straight forward. In order to complete the trials with the two-inch 
diameter LDPE, the 94 cm (37 in) section of the tow ta nk in front of the high speed 
camera was isolated with two wooden rods. The LDPE sections were placed into three 
different groups. Each group contained 54 rods, which equated to a 0.1667 ice density 
ratio for the tow tank segment. When all rods were added to the tow tank section, it was 
at 50% capacity. The term capacity is used because it is impossible to fully cover a 
rectangular area with circles. A 100% capacity would be the maximum number of LDPE 
rods that could be placed inside the ice section. The actual percent area of the section 
covered by the ice equivalent at 50% capacity is 37.5%. A rectangular pattern of circles 
was used for this calculation. All trials for the various hull orientations were completed 
first for the lower ice density. Time was allowed between trials for the tow tank water 
and LDPE rods to reach a state of rest. Upon the completion of the low ice density ratio 
testing, the additional LDPE segments required for higher ice density ratio were added to 

the tow tank. The experimental setup before each trial is shown in Eigure 38. 
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Figure 38. Setup for 50% Capaeity Small Ice Trials. 

Next, the four-inch diameter LDPE chunks were added to the tow tank in 
accordance with the test plan. In order to ensure the LDPE segment remained in place, 
holes were drilled near the top of the rods. A small loop of soldering wire was inserted in 
one of the rods for the single strike tests. Two other rods were connected with a small 
length of soldering wire at the desired separation length. The wire was loosely placed in 
the second rod so the bow wave would separate the two rods from each other. Soldering 
wire was used because it deforms easily. The attachments are shown in Eigure 39. 
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Figure 39. Soldering Wire Attachments to Aid in Stabilizing Large Ice 

Equivalent. 


A length of thread was attached to the soldering wire. The thread was then 
connected to the pulley system, which was located along the centerline of the tow tank. 
Once the rods came to a state of rest at the center of the tank, the line was removed from 
the pulley system. The dangling thread was assumed to have a negligible effect on the 
experiment. The rod had to be manually placed for the off-centered single strike impact 
test. If the footage showed the rod moved out of location, the trial was repeated. The 
setup for the large ice impact tests is shown in Figure 40. 
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Figure 40. Large-Diameter Impact Test Setup. 


Once all of the experimental data was obtained, a Go-Pro was attached to the 
bottom of the plate. The Go-Pro was used to record underwater footage of the area 
immediately preceding the plate. The goal was to capture the interaction between the 
bow wave of various hull types and the ice equivalent in the plate’s path. The 
composite plate with the Go-Pro is shown in Figure 41. 



Figure 41. Go-Pro Used to Record Underwater Footage. 
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III. RESULTS 


A. EXPERIMENTAL DATA REDUCTION 

This experiment resulted in 104 Excel spreadsheets of raw experimental data. The 
data for each trial consisted of the force required to pull the carriage and the values 
recorded by the three strain gauges. First, there was a permanent bias in the force sensor 
attached to the front of the carriage. This was removed by subtracting the bias from the 
data sets. The tension experienced by the load sensor was recorded as a negative value. 
For this reason, the force data was also inverted. An example of the raw for and force 
data is shown in Figure 42. 



Figure 42. Raw Experimental Force Data. 

This data shows the force is originally offset by 40 Newtons (N). The purpose of 
the base water trials was to measure the average force required to pull the plate through 
the water for different speeds and orientations. This was obtained by taking the average 
of the steady state section of the data defined in Figure 12. The corrected experimental 
data for the trial shown in Figure 42 is displayed in Figure 43. Smooth lines were used to 
make the data curves more visible. 
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Figure 43. Corrected Experimental Force Data Used in Analysis with Steady 

State Region Isolated. 

The strain data also had a permanent bias for each strain gauge that had to be 
accounted for. The same process used on the force data was used on the strain data. 
Typical corrected strain data is shown in Figure 44. Figure 45 shows the location of the 
strain gauges and the terminology used to describe the strain location being discussed. 



Figure 44. Corrected Experimental Strain Data Used in Analysis. 
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Figure 45. Strain Gauge Loeation and Terminology. 

B. BASE WATER TRIALS 

While not the focus of this study, the performance of the semi-submerged plates 
for the three orientations was investigated. These trials show that it required nearly half 
the amount of force to move the conventional plate through the tow tank than the force 
required to move the other two plates. The average force required to move each plate 
orientation is shown in Figure 46. 
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Figure 46. Average Foree Required to Tow Different Plate Orientations for 

Base Water Trials. 

The disparity in the towing forees most likely stems from the flow eharacteristics 
of the water around each plate orientation. When the tow tank was at rest, 12.7 cm (5 in) 
of the composite plate was above the waterline. The length of the plate above the 
waterline will be mentioned as the freeboard. The conventional oriented plate deflects the 
water down, whereas the tumblehome plate deflects the water up. This creates a large 
bow wave in front of the tumblehome hull. The water also rises higher on the 
tumblehome plate as it travels faster through the water resulting in little to no freeboard. 
As the plate velocity increased, the water partially topped the tumblehome oriented plate. 
This results in the tumblehome plate having to displace more water, increasing the force 
experienced by the load sensor attached to the pulley system. On the other hand, most of 
the five inches above the waterline remained dry as the conventional plate was towed 
through the water. The vertical plate exhibited the middle ground between the two 
extremes of the tumblehome and conventional flow profiles. The flow around the vertical 
plate is shown in Figure 47. This shows that the amount of freeboard was slightly reduced 
as the plate moved through the water. Figure 48 shows that water flowing over the 
tumblehome plate. Figure 49 shows the freeboard does not appreciably change for the 
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conventional plate. These images show the turbulent nature of the water flowing around 
the vertical and tumblehome oriented plates when compared to the flow around the 
conventional oriented plate. 



Figure 47. Flow Around Vertical Plate at 9 Hz Exhibiting Partial Freeboard. 



Figure 48. Flow Around Tumblehome Plate at 9 Hz Exhibiting Overflow. 
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Figure 49. Flow Around Conventional Plate at 9 Hz Exhibiting Full Freeboard. 

The flow profile around each plate orientation is also responsible for the strain 
experienced by the composite plate. The tumblehome orientation resulted in the greatest 
strain experienced by the plate, while the conventional orientation resulted in the lowest 
strain experienced by the plate. As predicted by the Ansys model, the center strain gauge 
experienced the most strain. This is a result of its close proximity to the fixed support. 
The strain recorded by the vertical strain gauges is displayed in Figure 50. 
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Figure 50. Average Strain Experienced by Vertical Strain Gauges during Base 

Water Trials. 


The strain experienced by the horizontal strain gauge was more complex when 
compared to the vertical strain. It is logical for the strain experienced by the plate to 
increase at it travels faster through the water. This is exhibited by the vertical strain 
gauges. The horizontal strain values reveal a different strain profile. This is shown in 
Figure 51. 
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Figure 51. Average Strain Experieneed by the Horizontal Strain Gauge for Base 

Water Trials. 

This strain profile is believed to be the result of two different modes of 
deformation eompeting with eaeh other. The two different modes of deformation are 
bending around the vertieal eenterline of the plate and bending around the horizontal 
eenterline of the plate. To begin, it is important to notiee that strain values in the 
horizontal direetion are between one and two orders of magnitude less than the strain 
experienced in by the vertical strain gauges. This is a product of the location of the fixed 
support and the waterline. Since most of the pressure against the plate is below the fixed 
support, the plate primarily experiences bending around the horizontal centerline. 
Bending along both centerlines gradually increases until the strain reaches a critical point 
near 0.9 m/s. At this point, the deformation along the horizontal centerline becomes the 
dominant mode. As deformation increases in one mode, it decreases the amount of 
deformation for the other. Figure 50 shows that the vertical strain is greatest for the 
tumblehome orientation. This means the plate is bending more around the horizontal 
centerline. This explains why the strain decreases the most in the horizontal direction for 
the tumblehome orientation at high velocities. Conversely, the vertical strain for the 
conventional plate orientation has the lowest strain. Therefore the strain in the horizontal 
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direction has the potential to be the highest value amongst the three orientations. The two 
modes of deformation are shown in Figure 52. The trade-off in deformation was also 
apparent in the small ice density tests where strain in the horizontal and vertical direction 
was inversely proportional. 



Vertical Bending Axis 

Horizontal Bending Axis 
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Figure 52. Competing Deformation Modes at High Velocities. 

C. SMALL ICE DENSITY TRIALS 

In order to assess how the presence of the two-inch diameter LDPE rods affected 
the strain experienced by the plate and the force required to pull the plate through the 
water, the collision event was isolated. The distance from the starting location of the 
carriage to the LDPE rods was measured. Using the velocities obtained from the high 
speed camera, the amount of time required for the composite plate to reach the section of 
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ice was calculated. This time was then overlaid on the experimental data for the strain 
and foree graphs. If the LDPE had an effect on the strain or foree, it should have been 
visible on the plotted graphs at, or immediately following the calculated time of impact. 
For example, using the earriage veloeity from the 7 Hz speed trial, it would take the 
carriage three seconds to reach the ice section. A three second window was created in the 
strain data. This window started at the first large strain spike at two seconds which 
indicates that the carriage began moving. The window stops just after five seconds. The 
spike in strain values that immediately follow the time window results from the 
interaction between the LDPE rods and the composite plate. An example of this proeess 
is shown in Figure 53. 



Figure 53. Else of Carriage Veloeity and Distance to lee Equivalent Seetion to 

Verify Impact Time. 

The section immediately after this window was then isolated in order to obtain the 
difference between the peak strain from the steady state strain experienced prior to the 
impact. The isolated section of strain data lies between the last solid vertical line and the 
dashed vertical line in Figure 53. This is shown in Figure 54. The peaks are less defined 
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than in those in Figure 53 beeause the time seale has been expanded nearly tenfold. This 
proeess was also eompleted for the foree values. 



It was diffieult to make eomparisons between the performanees of each hull 
orientation as they interact with the ice equivalent. The steady state strain and force 
values for each plate orientation are different. Since the values are different, the 
performance of each plate orientation was assessed by the percent increase or decrease 
between the steady state value (strain or force) prior to the ice impact and the peak value 
that resulted from the collision with the ice equivalent. Trials at 3 Hz were excluded from 
this analysis because the variation in the strain and force values was too inconsistent. 

1, Peak Force Analysis 

The analysis of the ice trials at 50% capacity yielded a common trend for the 
percent increase in the peak force required to tow the composite plate. The percent 
increase in force for the different plate orientations is shown in Figure 55. 
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Figure 55. Peak Force Percent Increase for High Capacity Ice Trials. 

As the velocity of the plate was increased, the percent increase in the peak force 
gradually decreased. This is due to the fact that the force required to move the carriage at 
high velocities is much higher than the force required to move the plate at low velocities. 
The force of drag on the plate is a function of velocity squared. While the percent 
increase in peak velocity diminishes, the increase in value is still notable. For instance, at 
a speed setting of 5 Hz (0.85 m/s), the average force required to pull the conventional 
oriented plate before it reached the ice section was 5.35 N. The 92% increase resulted in a 
peak force of 10.3 N, a rise of 4.95 N. On the other hand, at 9 Hz (1.51 m/s) the average 
force was 45 N. The 34% increase in force increased the value to 61 N, an increase of 
16 N. That increase is more than three times that of the lower speed trial. 

The data shows that the tumblehome orientation experienced the smallest increase 
in peak force values. This initial study shows that a tumblehome hull travels more 
efficiently through Arctic sea lanes congested with free-floating ice. The footage 
recorded by the high-speed camera supports the fact that the tumblehome orientation 
navigates better in the presence of ice. Still images of the footage were taken at three 
locations along the ice section in the tow tank. The first set of images shown in Figure 56 
display the three plate orientations as they approach the ice section. 


58 






















Figure 56. Profile View of the Small lee Seetion as the Three Plate Orientations 

Approaeh the lee Equivalent. 


While the flow patterns around the plate were not drastieally different to begin 
with, slight differenees are notieeable. For instance, as the tumblehome-oriented plate 
approached the ice equivalent, the ice equivalent was pushed away from the plate. Not 
only was it pushed away in the lateral direction, but the rods began to rise in the bow 
wave generated by the plate. The bulbous bow of large navy ships most likely creates a 
bow wave similar to that of the tumblehome plate, which displaced the ice equivalent. 
However, small craft such as the M-80 Stiletto are not large enough to constitute the 
addition of a bulbous bow. The ice equivalent did not experience much displacement for 
the other two plate orientations. The next sets of images in Figure 57 were taken as the 
plate contacted the ice equivalent. 


59 






Figure 57. Profile View of the Small lee Section as the Three Plate Orientations 

Contact the Ice Equivalent. 


Each image has a red arrow which shows the direction in which the ice equivalent 
was deflected once it contacted the plate. The image of the vertically oriented plate shows 
perturbations in the water in front of the plate highlighted by the black arrow. This is the 
result of the individual ice equivalent segments impacting one another. Unlike the other 
two plate orientations, the vertical plate orientation does not deflect the ice equivalent up 
or down. This causes the ice equivalent in the front of the section to be driven into the ice 
equivalent behind it. The conventional-oriented plate deflected the ice equivalent deeper 
into the water and underneath the plate. This may be problematic for sea-going vessels as 
the ice may deteriorate the quality of the propellers over time, resulting in a decrease in 
performance and an increased noise signature. On the other hand, the tumblehome- 
oriented plate pushed the ice equivalent up. The difference between these two flow 
profiles demonstrates why the tumblehome oriented plate experienced a smaller 
percentage increase in peak force when compared to the conventional oriented plate. 
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Rather than colliding with the ice equivalent, the tumblehome orientation displaces the 
ice around the plate. The last set of images in Figure 58 show various plate orientations 
near the end of the ice section. 






Conventional 


Tumblehome 


Figure 58. Profile View of the Small Ice Section as the Three Plate Orientations 

near the End of the Ice Section. 


These images provide further evidence of the findings mentioned above. The 
image of the conventional orientation shows a large amount of ice equivalent being 
displaced deep into the water. On the contrary, there appears to be a large concentration 
of ice equivalent remaining in front of the tumblehome oriented plate. The ice equivalent 
also remains near the water level for the vertical plate. It appears that nearly half of the 
ice equivalent has flowed around the vertical plate. Two enlarged images that display the 
flow profile of the ice equivalent around the tumblehome-oriented plate are shown in 
Figure 59 and Figure 60. These better show the flow of the ice equivalent around the 
plate. 
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Figure 59. Enlarged Image of the Flow Profile for Tumblehome Plate in the 

Middle of the lee Section at 9 Flz. 



Figure 60. Enlarged Image of the Flow Profile for the Tumblehome Plate at the 

End of the Ice Section at 7 Hz. 
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2, Peak Strain Analysis 

The composite plate experienced the same phenomena mentioned in the base 
water section where the bending around the horizontal and vertical axis was inversely 
related. The strain values recorded from each gauge was expected to increase when the 
plate collided with the ice section. Instead, the horizontal strain increased while the 
vertical strain values decreased. This behavior is attributed to the competing deformation 
modes explained in the base water trials section. The collision with the LDPE rods 
caused the plate to experience increased bending along the vertical centerline. This was 
because the ice equivalent was present at the waterline. When the bow wave was 
accounted for, the ice equivalent at the waterline added pressure against the plate near the 
top of the plate where the fixed support was present. This added pressure resulted in 
bending about the fixed support. The increased deformation around the vertical centerline 
decreased the amount of bending around the horizontal centerline. The peak strain in the 
horizontal direction always increased. The data is presented to allow for the comparison 
between the peak strain values experienced by the different plate orientations. Figure 61 
shows the percent change in strain experienced by the vertical rosette strain gauge. The 
220% increase in the plot was due to the low value of the average strain before the impact 
event. The strain recorded from the horizontal strain gauge increased from 50 ps to 
160 ps, while the vertical strain at the same location decreased from 160 ps to 100 ps. 
The average strain was in the denominator for the percent change calculation. 
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Figure 61. Percent Change in Strain for the Vertical Orientation during High- 

Capacity Small lee Trials. 

The 7 Hz (1.18 m/s) speed trial (center data point) for the vertical plate resulted in 
the largest percent change for all strain gauges. This trend is also present for the data 
obtained in the tumblehome trials. This is most likely due to the bow wave generated in 
front of the composite plate for the vertical and tumblehome plate orientations. The 
images of the different oriented plates as they approach the ice section shown in Figure 
56 display this bow wave. The plate velocity affects the interaction between the plate and 
ice equivalent in two ways. First, the force resulting from impact is directly related to the 
speed at which the plate strikes the ice. Increasing impact velocities result in greater 
incident forces. The velocity of the plate also affects the size of the bow wave generated. 
Increasing plate velocity results in an increased bow wave. This bow wave aids in 
displacing the ice equivalent, lowering the initial impact between the plate and ice 
equivalent. There is a critical point for the plate velocity where, at some speed slightly 
higher than the 7 Hz speed setting , the diminishing effect on the incident force from the 
bow wave generated by the plate surmounts the amplification in force resulting from the 
impact at higher speeds. The percent change in strain for the tumblehome orientation is 
shown in Figure 62. 
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Figure 62. Percent Change in Strain for the Tumblehome Orientation during 

High-Capacity Small Ice Trials. 

The dominant bow wave effect is apparent between the second and third data 
points for the three strain gauges. The strain experienced by the conventional oriented 
plate is unique. Since the bow wave generated is smaller, its diminishing effects appear at 
lower velocities. With the exception of one trial, the high-speed trial resulted in the 
greatest increase in strain. The percent increase in strain for the conventional oriented 
plate is displayed in Figure 63. 
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Figure 63. Percent Change in Strain for the Conventional Orientation during 

High Capacity Small Ice Trials. 

Some trials displayed a dynamic response after the initial impact. The resulting 
strains values appeared to be random in nature and no trends could be determined. The 
dynamic response can be explored in a future experiment in which the section of ice 
extends for a larger percentage of the tow tank. In this experiment, the ice section was 
only 94 cm (37 in) long. If the post-impact dynamic response is desired, the ice section 
should extend at least three meters. A three meter ice section at 50% capacity would 
require 531 two-inch diameter LDPE rods. While this would complicate running multiple 
trials, it would be necessary to fully characterize the performance of each hull orientation 
in the Arctic environment. 

3, Ice Capacity Effect 

This experiment showed that ice capacity does not have a large effect on the 
initial peak in strain values produced from the impact of the plate and the ice section. All 
trials were completed at 7 Hz (1.18 m/s). The largest percent change in strain for all trials 
was 2.2%. The largest percent change for the center strain gauge was 0.16%. Once again, 
the percent change in strain for the vertical and tumblehome orientations have similar 
trends, while the conventional orientation has a unique data trend. An increasing ice 
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capacity yields greater percent ehange in strain for the vertieal and tumblehome 
orientations. The pereent ehange in strain for the varying plate orientations is shown in 
Figure 64, Figure 65, and Figure 66. Most values of percent ehange were negligible. 



Figure 64. Pereent Change in Strain for the Vertical Orientation at 7 Hz for 

Different lee Capacities. 



Figure 65. Pereent Change in Strain for the Tumblehome Orientation at 7 Hz 

for Different Ice Capaeities. 
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Figure 66. Percent Change in Strain for the Conventional Orientation at 7 Hz 

for Different Ice Capacities. 


D, LARGE ICE TRIALS 

After viewing the data obtained from both the single impact and double impact 
large ice trials it was determined that, with one exception, a single impact event with the 
large diameter ice equivalent could only be isolated for the 7 Hz (1.18 m/s) speed trials. 
The pulley system used to translate the carriage across the tow tank creates too much 
variation in the recorded values of force and strain. The impact event is concealed in the 
fluctuating force and strain values for low carriage speeds. This is shown in Figure 67. 
The vertical line indicates the time at which the impact should have occurred. The strain 
recorded by the center strain gauge was added to the plot. The strain values aid in 
visually identifying the moment of impact. 
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Figure 67. 5 Hz Conventional Trial with no Apparent Impact Point. 


The fact that only the trials from the 7 Hz yielded noticeable results may not be by 
chance. As discussed for the small ice trials, 7 Hz was the critical speed at which the bow 
wave had yet to diminish the impact between the plate and ice equivalent. This also 
pertained to larger sized ice equivalent. The exception to this fact is the high speed 
conventional trial. This is due to the insignificant bow wave created by the conventional 
orientation. The difference between the data obtained from a 7 Hz trial and a 9 Hz trial is 
shown in Figure 68 and Figure 69, respectively. The spike in strain is most evident. 
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Figure 68. Strain and Force Resulting from Large Ice Double Impact at 7 Hz. 
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Figure 69. Strain and Force Resulting from Large Ice Double Impact at 9 Hz. 


The percent change in stress and strain for the 7 Hz trials is tabulated in Table 6. 
The spike in force was greatest for the conventional orientation. The bow wave was not 
large enough to begin displacing the ice equivalent. This would have softened the impact 
between the plate and ice equivalent. 
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Table 6. Percent Change for Force and Strain Values at 1.18 m/s Due to Impact 

with Large Ice Equivalent. 


Orientation 

Force 

Vertical Rosette 

Horizontal 

Rosette 

Vertical Center 

Vertical 

6.68% 

-25.00% 

140.00% 

-11.84% 

Conventional 

42.70% 

-41.67% 

85.71% 

-16.95% 

Tumblehome 

15.01% 

16.00% 

100.00% 

-8.60% 


The results in Table 6 are best explained by inspecting the still images taken 
during the large ice trials. These images show that the relationship between the different 
plate orientations and the resulting behavior of the ice equivalent are similar to those of 
the small ice trials. The only difference between the two is that the larger rods have more 
inertia. This increase in inertia lessens the pushing effect of the bow wave preceding the 
plate. The images are presented just as they were presented in the small ice section. The 
first set of images in Figure 70 shows the composite plate as it first approaches the large 
ice equivalent. The plate in the tumblehome orientation displaces the ice equivalent 
upwards and away from the plate. Conversely, the ice equivalent experiences little to no 
displacement for the conventional plate orientation. The next set of images in Figure 71 
show the profile view immediately following the impact between the plate and the ice 
equivalent. The conventional oriented plate deflects the ice equivalent deeper into the 
water, whereas the ice equivalent moves with the tumblehome oriented plate. The ice 
equivalent also remains in front of the vertical oriented plate. The final set of images in 
Figure 72 show the continuation of each process. The ice equivalent continued to travel 
in front of the plate for the vertical and tumblehome oriented plates. The interaction 
between the ice equivalent and the plate in these orientations is prolonged. This results in 
a low force impact for a long length of time. On the other hand, the conventional 
orientation results in a nearly instant, high force impact. 
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Figure 70. Profile View of each Plate Orientation on Approach to 

the Large Ice Equivalent. 



Figure 71. Profile View of each Plate Orientation on Contact with 

the Large Ice Equivalent. 
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Figure 72. Profile View of Eaeh Plate Orientation After Contact with 

the Large Ice Equivalent. 


E. FRAMED COMPOSITE PLATE 

For the plate fixed to the attachment bar, there was a notable difference between 
the strain at the center of the plate and at the offset location. This difference was not as 
evident when the frame was used to attach the plate to the attachment bar. The average 
strain in the vertical direction is shown in Figure 73. The maximum value in strain was 
still located at the center of the plate. The data from the center strain gauge is designated 
by the dashed lines. 
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Figure 73. Average Strain Experienced by the Framed Plate during 

Base Water Trials. 

As predicted, the change from the fixed support of the attachment bar to the fixed 
boundary conditions from attaching a frame to the composite plate greatly reduced the 
strain experienced by the plate. The strain values ranged from 13 ps to 1280 ps for the 
plate attached to the attachment bar. The framed boundary condition resulted in this range 
to be decreased to 14 ps to 325 ps. The same pattern occurred for the horizontal strain 
gauge. 

The interaction between the ice equivalent and the framed plate was visible in the 
strain data. For the most part, there was a response in the strain gauges at the time of 
impact. However, it would not make sense to compare the percent increase or decrease in 
strain to the frameless plate. Since the strain values are very small, a small change in 
strain would appear to have a large impact on the percent change in strain. For example, 
for the framed plate at 7 Hz, a change of 10 ps would result in a 10% change in strain 
value. If the strain of the frameless plate was changed by 10 ps, the resulting percent 
difference would only be 2.4%. The difference in the location of the fixed support results 
in the two plates having different modes of deformation. For this reason, it makes little 
sense to compare the two. 
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F. 


SUMMARY OF EXPERIMENTAL RESULTS 


The results of this experiment were highly dependent on the location of the 
waterline of the towed composite plate in relation to the location of the fixed support 
pulling the plate through the tank. The composite plate was originally given 12.7 cm 
(five inches) of freeboard in hopes that the water would not flow over the top of the plate 
for high speed trials. With the exception of the 9 Hz (1.51 m/s) tumblehome trials, the 
water did not flow over the plate. The location of the waterline was important because the 
goal of this thesis was to evaluate the performance of different hull orientations as they 
travel through ice-congested waterways. The interaction between the plate and the ice 
occurs at the surface of the water. Figure 74 shows the spatial orientation of the fixed 
support of the plate and the pressure generated by the impact between the plate and ice 
equivalent. The yellow area signifies the ice equivalent along the waterline. 


I I Fixed Support 


n 



Figure 74. Incident Pressure Resulting from the Plate Travelling 

through the Ice Equivalent. 

The resulting force of the ice equivalent against the plate will incur bending about 
the vertical centerline. After observing the tow tank trials, it was apparent that the bow 
wave increased the height of the location of the ice equivalent as it impacted the plate. A 
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more realistic impact plane is designated by the black dashed line. The competing 
deformation modes of the plate were discussed heavily in the results section. Any 
deformation in the horizontal direction resulted in decreased deformation in the vertical 
direction and vice versa. This was first evident in the base water trials where increased 
vertical deformation decreased the amount of deformation recorded by the horizontal 
strain gauges. The same phenomena occurred during ice density trials. The presence of 
the ice equivalent resulted in increased bending about the vertical centerline, reducing the 
strain recorded by the vertical strain gauges. The results of this experiment would have 
been very different had the fixed support been place higher on the plate. This is portrayed 
in Figure 75. 



Figure 75. Incident Pressure against the Plate for the Case 
where the Water Line is Lowered. 

If this was the case, the plate would experience little to no strain in the horizontal 
direction. The impact force of the ice equivalent would further increase the strain in the 
vertical direction. All of this proves the dependability of the strain experienced by the 
plate on the physical location of the waterline on the plate and the location of the fixed 
support. 
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The strain data is unique to this experiment. For this reason, the most valuable 
information obtained during this experiment was the effect of the ice on the force 
required to tow the composite plate and the footage obtained that show the interaction 
between each plate orientation and the ice equivalent. Both of these showed that the 
conventional plate orientation performed the worst in the presence of ice. The still images 
of the footage revealed a lack of bow wave generated by this orientation. This resulted in 
an immediate and violent impact between the plate and the ice equivalent. The 
tumblehome oriented plate experienced the lowest percent increase in peak force incurred 
by the impact. The video demonstrated that the ice equivalent was displaced up and 
around the plate due to the bow wave, resulting in little to no impact. This bow wave may 
still be present for large navy ships with conventional hulls due to their bulbous bow. 
However, small craft operating in Arctic sea lanes are not large enough for bulbous bows 
and will experience flow patterns similar to this experiment. Furthermore the footage of 
the conventional oriented plate revealed that after the impact, the ice equivalent is forced 
underneath the plate, pushing it deep into the water. This could result in the ice being 
ingested by the propellers. If this occurs for a prolonged amount of time, the propellers 
could experience a decrease in performance and an increase in the noise signature of the 
vessel. One potential problem with the flow profile of the ice around the tumblehome 
oriented plate is the accumulation of ice in front of the plate. This build up was apparent 
for the small ice equivalent section used in this experiment. This is shown in Figure 76. It 
is possible that if enough ice is accumulated in front of the plate it could impede the 
motion of the plate. If the accumulation was severe, the ice would grind against the 
composite. This grinding motion could result in the failure of the fibers of the composite. 
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Figure 76. Ice Accumulation Present Preceding the 
Tumblehome-Oriented Plate. 


Another outcome of this experiment was the validation of LDPE as an ice 
equivalent. Throughout testing, the plastic maintained the same density and did not 
absorb any water. It allowed for repeatable experimentation and was simple to maintain 
and use. This method forgoes the extensive energy requirement of ice tanks that maintain 
the working area at freezing temperatures. 

1, Recommendations for Future Work 

If further information is sought regarding the flow of ice around moving 
structures, this methodology is recommended. Special consideration should be taken 
while considering the boundary conditions of the structure being towed through the tow 
tank. Testing on a prototype composite hull would be ideal. If another plate is used, a thin 
frame boundary condition would provide the most insight since most hulls are 
constructed of plates connected to the framework of the ship. A new frame should be 
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constructed, as the frame used in this experiment was too bulky and redueed the area of 
the plate exposed to the water. 

Future experimentation should be eonducted to better charaeterize the dynamie 
response of the plate after the initial impaet with the iee seetion. To aeeomplish this, the 
ice seetion must be extended mueh farther than the one in this experiment. A three-meter 
iee seetion should allow for a full analysis of the response following the initial impaet. 
A prolonged ice seetion would also answer the question of whether or not the iee 
aeeumulation effeet exhibited by the tumblehome orientation would negatively influence 
its performanee. It may also lead to further insight into the performanee of the other plate 
orientations. This testing should be eondueted with both small and large diameter LDPE 
rods. This experiment showed that the inertia of the iee equivalent has an effeet on the 
interaetion between the iee equivalent and the plate. If a longer iee seetion is used, an 
efficient method used to eolleet the ice equivalent would be necessary. The experiment 
ean also be repeated using higher ice eapaeity values to simulate more eongested sea 
lanes. 
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IV. CFX MODEL 


A CFX model of the experiment was created in an attempt to replicate the 
experimental data. Ideally, the CFX model would have been coupled to a structural 
model. This would allow for the pressure against the test plate calculated by the CFX 
model to be imported into a structural analysis tree. In order to do this there must be 
nodal connectivity between the fluid and the structure. This means the mesh of the fluid 
domain must be the same as the mesh of the structural analysis. This analysis has been 
completed for a fully submersed structure [26]. This CFX analysis was more complicated 
as it involved a semi-submerged plate. This CFX simulation required a mesh refinement 
in order to accurately model the free surface effect along the air-water interface. Without 
the mesh refinement, the boundary between the air and water would become indistinct. In 
order to assess the performance of each orientation, the incident fluid pressure on the 
plate was determined. All CFX results were steady state solutions. In the end, the 
presence of ice could not properly be modelled in the CFX analysis. The model presented 
does accurately show the free surface effect of the two-fluid domain around the plate for 
the conventional and tumblehome orientations. Ansys tutorials [27] were used as the 
basis for creating the two-fluid domain. This analysis was done for the conventional and 
tumblehome hull types. 

A, CONTROL VOLUME 

The control volume used for the CFX simulation was created in Solidworks. The 
control volume had the same width and height of the tow tank, but was shortened in order 
to keep the amount of mesh elements relatively low. Even with the reduction of the 
meshing elements, each CFX run required around 13 hours to meet the solver control 
requirements defined later. After trial and error it was determined that the control volume 
would extend 1.52 m (five feet) upstream of the plate and 3.048 m (10 feet) downstream 
of the plate for the angled orientations. This volume allowed for ample space to visualize 
the flow characteristics around the plate. The control volume for the vertical plate was 
3 m (9.84 feet) long. 
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The test plate was ereated by creating a void inside of a rectangular prism (the 
tow tank volume). The first step in creating the control volume was to create a rectangle 
with the dimensions of the tow tank. The plate orientation was then established by 
creating a plane at the desired plate location. The plate orientation is shown in Figure 93 
located in Appendix C, which contains all of the supplemental information regarding the 
Solidworks and CFX processes. The plate was then drawn and extruded from this plane. 
The “Merge Results” box was checked when extruding the attachment bar from the plate. 
The attachment bar was extruded a distance that allowed it to protrude from the 
rectangular prism of the tow ta nk once the rectangle was extruded. The rectangle was 
then extruded 4.572 m to achieve the necessary control volume in front and behind 
the test plate. Finally, the plate and attachment bar were removed from the control 
volume by selecting ‘Tnserf’ from the drop down menu and selecting “Features...”, 
“Combine.” The subtract feature was selected and the plate was removed from the control 
volume. This is shown in Figure 93 of Appendix C as well. The control volume was then 
imported into CFX. This process was completed for the three different hull orientations 
used in the experiment. The control volume for the tumblehome oriented plate is shown 
in Figure 77. 



Figure 77. Control Volume with Location of the Inlet, Outlet 

and Opening Displayed. 
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The three eontrol volumes were then meshed using the sizing parameters shown 
in Figure 78. The resulting mesh parameters of each control volume are shown in Table 
7. An image of the mesh used can be found in Figure 94 of Appendix C. Tetrahedron 
mesh elements were used. 


Sizing 

Use Advanced Size Function 

On: Curvature 

Relevance Center 

Coarse 

Initial Size Seed 

Active Assembly 

Smoothing 

Medium 

Transition 

Slow 

Span Angle Center 

Fine 

Cunrature Normal Angle 

3.0* 

Min Size 

2.5e.002 m 

Max Face Size 

2.5e.002 m 

Max Size 

2.5e.002 m 

Growth Rate 

Default (1.20) 

Minimum Edge Length 

2.e4)02 m 

L 

Use Automatic Inflation 

None 

Inflation Option 

Smooth Transition 

Transition Ratio 

0.77 

Maximum Layers 

5 

Growth Rate 

1.2 

Inflation Algorithm 

Pre 

View Advanced Options 

No 


Figure 78. Mesh Sizing Used for CFX Control Volume. 


Table 7. Mesh Parameters of Control Volumes. 


Plate Orientation 

Vertical 

Conventional 

Tumblehome 

Nodes 

237,599 

611,014 

610,977 

Elements 

1,361,902 

3,535,883 

3,535,750 


1, Fluid Parameters and Boundary Conditions 

Next, the boundary conditions of the fluid domain were defined. First, variables 

were created that were used to define fluid properties such as water height, densities, 

pressures, and volume fractions between the air and water. The full list of expressions 

created and their values is shown in Figure 95 in Appendix C. The water heights were 

fine-tuned in order to achieve the proper water level against the plate identical to that of 

the experiment. All fluid velocities were set to one m/s. 
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The two-fluid domain was created under the default domain tree, which highlights 
the entire control volume. Both air (at 25°C) and water were defined. Buoyancy was 
accounted for in the model in order for the air and water to separate into two distinct 
bodies. The Basic Settings of the fluid domain is found in Figure 96 in Appendix C. 
Table 8 shows the fluid models used. The boundary conditions of the control volume are 
defined in Table 9. Table 10. Table 11. and Table 12. 


Table 8. Fluid Domain Settings Establishing Two-Fluid System, after [27]. 


Tab 

Setting 

Input 

Fluid Models 

Multiphase 

Homogenous Model 

Checked 

Free Surface Model 

Standard 

Fluid Models 

Heat Transfer 

Homogeneous Model 

Checked 

Fluid Temperature 

25 [Cl 

Fluid Models 

Turbulence 

Option 

k-Epsilon 

Wall Function 

Scalable 

Fluid Pair Models 

Air Water 

Surface Tension Coefficient 

Checked 

Surf. Tension Coeff. 

0.072 [Nm^-l] 

Fluid Pair Models 

Surface Tension Model 

Option 

Continuum Surface Eorce 

Primary Fluid 

Water 

Fluid Pair Models 
Interphase Transfer 

Option 

Eree Surface 

Mass Transfer Option 

None 

Initialization 

Static Pressure 

Option 

Automatic With Value 

Relative Pressure 

UpPres 

Initialization 

Turbulence 

Option 

Medium (Intensity = 5%) 

Initialization 

Volume Fraction-Air 

Option 

Automatic with Value 

Volume Fraction 

UpVEAir 

Initialization 

Volume Fraction-Water 

Option 

Automatic with Value 

Volume Fraction 

UpVEWater 
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Table 9. Inlet Boundary Conditions, after [27]. 


Tab 

Setting 

Input 

Basie Settings 

Boundary Type 

Inlet 

Loeation 

Inlet 

Boundary 

Details 

Normal Speed 

1 [ms^-1] 

Turbulenee 

Intensity and Length Seale 

Fraetional Intensity 

0.05 

Eddy Length Seale 

UpH 

Fluid Values 

Air 

Value-Up VL Air 

Water 

Value-UpVLWater 


Table 10. Outlet Boundary Conditions. 


Tab 

Setting 

Input 

Basie Settings 

Boundary Type 

Outlet 

Loeation 

Outlet 

Boundary 

Details 

Llow Regime 

Subsonie 

Mass and Momentum 

Option-Average Statie 
Pressure 

Relative Pressure 

DownPres 

Pres. Profile Blend 

0.05 


Table 11. Opening Boundary Conditions, after [27]. 


Tab 

Setting 

Input 

Basie Settings 

Boundary Type 

Opening 

Loeation 

Opening 

Boundary 

Details 

Mass and Momentum 

Option-Entrainment 

Mass and Momentum 

Relative Pressure - 0 [Pa] 

Turbulenee - Option 

Zero Gradient 

Lluid Values 

Air 

Value-1 

Water 

Value-0 
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Table 12. Wall Boundaries. 


Tab 

Setting 

Input 

Basie Settings 

Boundary Type 

Wall 

Loeation 

Walls 

Boundary 

Details 

Mass and Momentum 

Option-Free Slip Wall 

Mass and Momentum 

Wall Adhesion-None 


2. Solver Control and Mesh Adaptation 

In order for CFX to eomplete the simulation, upwind differeneing had to be used. 
The other options resulted in the error messages part-way through the simulation. A mesh 
adaptation was added to better define the interfaee between air and water. This was 
required sinee adding a semi-submerged plate into the flow greatly inereased the 
eomplexity of the flow field. Later images show the high resolution of this boundary. The 
solver eontrols are defined in Table 13. The added mesh refinement ean be found Figure 
97 in Appendix C. 


Table 13. Solver Controls, after [27]. 


Tab 

Setting 

Input 

Basie Settings 

Adveetion Seheme 

Option-Upwind 

Turbulenee Numeries 

Option-First Order 

Convergenee Criteria 

Type-RMS 

le-06 

Advaneed Options 

Dynamie Model Control 

Cheeked 

Multiphase Control 

Surfaee Tension Coeffieient 

Volume Fraetion Coupling 

Option-Coupled 


B, CFX RESULTS AND DISCUSSION 

All of the simulations were eompleted at 1 m/s. For this reason, the flow profile 
around the plate is tamer in the simulations when eompared to the flow profiles visible in 
the images taken during the experiment. The waterline in the two-fluid domain was 
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highlighted using an isosurface. The isosurface was placed at the “Water.Volume 
Fraction” value of 0.99. A vertical plane was also created that ran the length of the 
control volume. This shows that the water was below the isosurface. Figure 79 shows the 
water level for the entire control volume for both orientations. 



liume Fraction 


I 


1+000 
7.50e-001 
5.00e-001 
2.50e-001 



-1.00e-015 


Iz.SOOe-OOl 

5.000e-001 

2.500e-001 

1.000e-015 


Fraction 



Figure 79. Water Level for Tumblehome (Top) and Conventional (Bottom) 

Flow Simulations. 


Next, the region before the plate was expanded to view the bow wave created by 
each orientation. As expected, the bow wave was larger for the tumblehome trial. It 
appears that there is a small region of air trapped between the plate and the water for the 
conventional trial. The wake also extends farther back for the conventional trial. The 
tumblehome oriented plate is fully submerged. The region in front of the plate is found in 
Figure 80. The wake region is shown in Figure 81. These images prove that the setup 
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used for the simulation was eorrect. For eomparison, the tumblehome trial ean be found 
on the left, while the eonventional trial ean be found on the right. 



Figure 80. Exploded View of Water Level Preeeding the Tumblehome (Left) 

and Conventional (Right) Oriented Plates. 



Ligure 81. Exploded View of the Wake Region behind the Tumblehome (Left) 

and Conventional (Right) Oriented Plates. 


The veloeity profile of the water was also inspeeted. The veloeity profile 
presented in Ligure 82 shows high veloeity regions to the left and right of the plate. The 
highest water veloeity this region was 1.272 m/s, which is 27% faster than the inlet flow. 
There is a low velocity region immediately behind the plate. This image shows water 
velocity vectors above the fluid domain as well. Ligure 79 shows that the lowest water 
volume fraction in the control volume is LOT O'* This would be the equivalent of the 
amount of water present in the air due to humidity. The presence of this extra fluid is 
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assumed to have negligible effeets on the results. Water veloeities were similar for both 
orientations. 



Figure 82. Water Veloeity in the Control Volume. Inlet Veloeities are 1.0 m/s. 


There were two desired results from this analysis. Sinee iee eould not be added to 
the system, the first goal of this initial Ansys analysis was to loeate the high pressure 
fluid in front of the plate that aided in displaeing the iee equivalent in the experiment. 
The seeond was to determine the ineident foree found on the plate. 

In order to examine the pressure in the fluid, horizontal planes were ereated in the 
eontrol volume at various heights along the plate. Horizontal planes had to be used 
beeause the hydrostatie pressure inereases as the depth inereases for water. Any 
eomparison using vertieal planes showed high pressures at the bottom of the tank. The 
next set of images eompares the pressure surrounding the plate. As before, the 
tumblehome trials ean be found on the left while the eonventional trials are on the right. 
Figure 83 shows the pressure ranging from the plane at the bottom to the plane at 
the surfaee of the water. The distanee between the heights of eaeh plane is 7 em 
(2.75 inehes). 
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Figure 83. Pressure Planes Signifying the High Pressure Region Preceding the 
Tumblehome (Left) and Conventional (Right) Oriented Plates. 
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The pressure preceding the tumblehome oriented plate was consistently greater 
than the pressure in front of the conventional oriented plate. In order to properly assess 
the pressure rise created by the plate, the hydrostatic pressure was subtracted from the 
pressure values from each plane. Table 14 shows the process of correcting the peak 
pressure value for each plane. It was assumed that the density of the water was 
1000 kg/m and that the highest plane was located at the water level. The latter 
assumption was made because the bow wave made it difficult to find the exact water 
level of the water. Regardless, the tumblehome trial consistently had greater pressure 
values than the conventional trial. The information was plotted in Figure 84. 


Table 14. Elimination of Hydrostatic Force from Peak Pressure Data. 


Plane 

Depth 

(cm) 

Resulting 
Hydrostatic 
Pressure (Pa) 

Pressure from 
Tumblehome 
Trial (Pa) 

Pressure from 
Conventional 
Trial (Pa) 

Final 

Tumblehome 
Pressure (Pa) 

Final 

Conventional 
Pressure (Pa) 

21 

2060.1 

2224.1 

2167.7 

164.0 

107.6 

14 

1376 

1936.3 

1654.6 

560.3 

278.6 

7 

686.7 

1205.6 

1094.7 

518.9 

408.0 

0 

0 

518.9 

454.5 

518.9 

454.5 



Figure 84. Pressure Rise Preceding Both Plate Orientations. 
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This increase in pressure may be the souree of the iee displaeing property 
exhibited by the tumblehome oriented plate in the experiment. The greatest pressure 
differential between the pressure in front of each orientation was near the lower section of 
the plate. 

Finally, the ineident pressure against the plate was determined. The data from the 
base water trials showed that it required nearly twiee the foree to tow the tumblehome 
oriented plate through the water than to tow the conventional oriented plate. The Ansys 
trials support this finding. The results of the Ansys trials shown in Figure 85 show that 
the incident pressure is consistently greater on the tumblehome oriented plate. The 
ineident force against the plate is the integral of the pressure over the surface area of 
the plate. 



[Pa] 


Figure 85. 


Ineident Pressure against the Tumblehome (Left) and Conventional 
(Right) Oriented Plate Aequired from Ansys Simulation. 


C. SUMMARY OF CFX ANALYSIS 

The CFX simulation properly eonstructed the tow tank environment. The two- 
fluid domain was established with the free surfaee present at the same location of the 
water line on the plate in the experiment. The high pressure region preeeding the plate for 
both the tumblehome and conventional orientation was successfully simulated. The 
pressure in front of the plate was greater at all water levels for the tumblehome-oriented 
plate. This increased pressure aided in the plate’s ability to displace the iee equivalent 
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during experimental trials. By displaeing the iee, the impaet between the plate and the iee 
equivalent was lessened. The ineident pressure against the plate was also plotted. Onee 
again, the simulation supported the findings in the experiment. In the experiment, the 
tumblehome oriented plate required a greater towing foree than the eonventional oriented 
plate required. The Ansys analysis determined that the pressure on the tumblehome 
oriented plate was greater than the pressure ineident to the eonventional oriented plate. 

1, Recommendations for Future Work 

While the high pressure region was suoeessfully simulated, adding the iee 
equivalent to the Ansys analysis was not aeeomplished. A eomplete analysis would be 
eompleted if two tasks were aeeomplished. First, future work ean be done to insert the iee 
equivalent into the Ansys simulation. It may be possible that Ansys is not eapable of 
aoeounting for solid free-floating partieles. If this is the ease, different modeling software 
suitable for FSI FEM must be used. Finally, in a eomplete analysis, a struetural analysis 
of a eomposite plate would be eoupled to the CFX analysis for the two tank experiment. 
This would allow for the ineident pressure on the plate to be imported to a struetural 
analysis to determine the stress and strain experieneed by the eomposite plate. 
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APPENDIX A. TOW TANK SOP 


A. FILLING AND DRAINING 

In order to operate the tow ta nk it is neeessary to understand the filling and 
draining proeure. The tow tank is shown below in Figure 86. 



Fill Valve Figure 86. Flalligan Hall Tow Ta nk . Drain Valve 


It is 36 feet long, three feet wide, and four feet high. These dimensions allow the 
tank to effeetively house a total of 2423 gallons of water. This ealeulation used a total 
water height of three feet rather than the ta nk height of four feet to allow for the waves 
ineurred during experimentation. In order to fill the ta nk the first step is to elose the drain 
valve. As shown in the figure, the drain valve is loeated on the right side of the tank. The 
open and elosed position of the valve is shown in Figure 87. 



Figure 87. Drain Valve Open (Left) and Closed (Right) Positions. 
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Once this valve is closed, the ta nk may be filled with water. The water line can be 
found on the opposite side of the tank from the drain valve. The open and elose positions 
are identieal to the drain valve. Copper piping has been added to the original system due 
to issues involving back flow, which resulted in the flooding of the laboratory floor. 
Rather than the water being added from the bottom of the tank where the water pressure 
is notably high once full, the water is redirected and is added from above the tank as 
shown in Figure 88. 



Figure 88. Tow Tank Fill Piping. 


B. CARRIAGE OPERATION 

With the ta nk ready for testing, the next step any experiment is to understand how 
the carriage operates. The carriage translates along the length of the tow ta nk along rails 
shown in Figure 89. The areas along the rails should be cleared of any objects. For safety 
reasons, when the tow ta nk is in operation, all personnel must stay behind the yellow 
floor tape surrounding the tank. A pulley system is used in eonneetion with the motors to 
translate the earriage along the tank. 
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Figure 89. Carriage System. 


Operators should exercise caution when the carriage is at the outer limits of the 
tow tank. There is no safety mechanism to stop the carriage from hitting the wall. This 
will cause major damage to the pulley system and its connections. First, the motor must 
be powered on using the outlet boxes located between the tow ta nk entrance of Halligan 
Hall and the drain valve. The motor switches are shown in Figure 90. The on position is 
clearly marked on the box. For powering on, the middle switch should be powered on 
before powering on the right box. For powering off the carriage system, simply reverse 
the order. 










The carriage is controlled with the controller shown in Figure 91. The controller 
has 100 feet of cable allowing its use along the entire tank. 



Figure 91. Carriage Controller. 


The important keys for operation are explained below: 

ENTER: Used to change speed settings. Different speeds are represented in Hz. Change 
the speed setting by using the up and down arrows, then pressing Enter again. 

LOCAL/REMOTE: Must be toggled until “Local” is displayed in the upper right corner 
of the controller. 

FWD/REV: Moves the carriage forward or backwards at the speed inputted. As of this 
thesis, REV moves the carriage from left to right. Carriage continues to move until 
Stop is pressed, 

JOG: Pressing Jog will illuminate a green light on the comer of the button. This feature 
allows the operator to move the carriage forward and backwards only when the 
EWD/REV button is depressed. The carriage will move at the 3 Hz speed setting only 
while the EWD/REV button is depressed under Jog operations. 

A summary for operation is found on the next page. A table with the resulting RPM of 
each speed setting is shown after the summary. 
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Operational Use-Summary 


Read Entire Page Before Operation 

• Ensure that the rails that the carriage moves on is clear of any objects 

• Power on the tank using the switches located near the entrance 

o Powering On - First flip middle switch, then outside switch (Labeled 
“Tow Tank”) 

o Powering Down - First flip outside switch, then middle switch 

• Retrieve the controller 

o Press “Local/Remote” until “Local” is displayed in upper right hand 
corner. 

o Press Enter to select a speed. Toggle up or down then press Enter again. 

o Ensure all personnel are outside the area cordoned off by the yellow 
caution tape on the floor. 

o Press Forward or Reverse to move carriage at selected speed. Carriage 
continues to move until Stop is pressed, 

o Another movement option is Jog. Jog moves the carriage only when the 
Forward or Reverse button is depressed. Press Jog (Green corner light will 
light up), then hold the desired direction for the carriage to move. 

• Power off the switches when finished with the tow tank 


Exercise caution near the outer limits of the tank. There are no safety 
mechanisms to stop the carriage from moving. A collision with the wall 
may break the pulley system or its connections. 
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Carriage Speeds in Water 


Entered Speed (Hz) 

Resulting RPM 

Carriage Speed (m/s) 

2.0 

N/A 

N/A 

2.5 

N/A 

N/A 

3.0 

74 

0.4995 

3.5 

88.5 

0.58415 

4.0 

103 

0.6688 

4.5 

117.5 

0.7534 

5.0 

132 

0.8381 

5.5 

146.5 

0.9228 

6.0 

161 

1.0074 

6.5 

175.5 

1.0921 

7.0 

188 

1.1767 

7.5 

200 

1.2614 

8.0 

214.5 

1.3460 

8.5 

227 

1.4307 

9.0 

240 

1.5153 
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APPENDIX B. STRAIN GAUGE CONNECTION: 
QUARTER-INCH BRIDGE ADAPTER 


Each time a new plate is tested in the tow tank, the strain gauges must be 
eonneeted to the data aequisition system. This seetion summarizes how to install a new 
set of strain gauges. The first step is to use a small flathead serewdriver to open the ports 
of eaeh quarter-ineh bridge adapter. The two wires from the strain gauges should be 
inserted to ports 0 and 2. The serews should be tightened to secure the wire inside the 
port. Finally, a bridge needs to be fabrieated from a small length of wire. This experiment 
used a one-ineh section of wire with its ends stripped. This small wire should eonnect 
ports 1 and 2. This wire ereates a ground for the voltage to travel to. To ensure the strain 
gauges are set up properly a voltmeter ean be used to eheek the resistanees aeross ports 
0 and 1. As long as there is a resistanee between the two ports it will funetion properly. 
The resistanee may not be the published resistanee of the strain gauge. Figure 92 shows 
the connections for a plate with three strain gauges. 



Figure 92. Strain Gauge Conneetion to the Quarter-Ineh Bridges. 
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APPENDIX C. SUPPLEMENTAL CFX IMAGES 


I Features | Sketch j Evaluate | Dir^pert^^^ Products 




S-® TankWPIateSweptSmallerT... 



^ X 


Operation Type 
Add 

o Subtract 




dl 0 C3 ^ 



%L»z 


Figure 93. The Towed Plate and the Attaehment Bar are Removed from the 
Control Volume of the Tow Ta nk Using the Combine Feature in 
Solidworks. The Plane Created to Generate the Plate Orientation is also 

Visible. 



ANSYS 

R15.0 


Figure 94. Mesh Used in CFX Trials. 
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d Expressions 



DenH 

(Den Water-DenRef) 

u 

DenRef 

1.185 {kg m^-3] 


DenWater 

997{kg m''-3] 


DownH 

0.73152{m] 


DownPres 

DenH^g'Vovi’n VFWater *(powr^-y) 


DownVFAir 

step((y-DowrH)/l {m]) 


DownVFWater 

1-DownVFAtr 


UpH 

0.73152{m] 


UpPres 

DenH^*UpVFWater*(UpH-y) 


UpVFAir 

step((y-UpH}/l {m]) 


UpVFWater 

1-UpVFAff- 


Figure 95. Variables Created for CFX Analysis, after [27]. 
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stalls of Default Domain in Flow Analysis 1 


Basic Settings 


Fluid Models | Fluid Specific Models | fluid Pair Models j Initialization | 


Location and Type 
Location 
Domain Type 
Coordinate Frame 


B158 

[Fluid Dom^ 

Coord 0 


□ 


Fluid ar»dPartide Definitions... B 


(E 


Water 0 



Domairt Models 



Air 

Water 



Figure 96. Fluid Domain Basic Settings, after [27]. 
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Dt<iiEoFM«b RerwiniK^t n Flow AKulvHf 1 
Sw? Stttpigs ; Adv»£«d OflbiKe. 

J AdlVfltS Adjpfton 
RfiQcnLijI: BlSS 

Save In wmeAaie Fls 
Adapben Cntcfia 


UbnotllkfUi^ 

tHr.VBltmt PractMin ^ 

Hu. 14 a. Steps 

2 

Optnn 

lUhplepflAbilHKh 

HiOtetCxir 

4 

Adapbon HrtFiod 


Op ton 

Sedition VbniiDtoi -r 


nnnvn tdoe Lengifv D.D 
Adaptnn Convcf gence CnlBia 


MaiK. 1 w. pet Step 


ReuiuaJTfpe 

RMS 

TvgclRcpdud 

D.Dai 


tocSciCi^ig Advanc)Hli<^lKn* _ 

Node Ak9c ParancEer 
hUiteiFfLcvcb i 

Figure 97. Mesh Adaptation Settings, after [27] 
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